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Abstract — The output of renewable energy fluctuates signif- 
icantly depending on weather conditions. We develop a unit 
commitment model to analyze requirements of the forecast output 
and its error for renewable energies. Our model obtains the 
time series for the operational state of thermal power plants 
that would maximize the profits of an electric power utility by 
taking into account both the forecast of output its error for 
renewable energies and the demand response of consumers. We 
consider a power system consisting of thermal power plants, 
photovoltaic systems (PV), and wind farms and analyze the effect 
of the forecast error on the operation cost and reserves. We 
confirm that the operation cost was increases with the forecast 
error. The effect of a sudden decrease in wind power is also 
analyzed. More thermal power plants need to be operated to 
generate power to absorb this sudden decrease in wind power. 
The increase in the number of operating thermal power plants 
within a short period does not affect the total operation cost 
significantly; however the substitution of thermal power plants 
by wind farms or PV systems is not expected to be very high. 
Finally, the effects of the demand response in the case of a 
sudden decrease in wind power are analyzed. We confirm that 
the number of operating thermal power plants is reduced by the 
demand response. A power utility has to continue thermal power 
plants for ensuring supply-demand balance; some of these plants 
can be decommissioned after installing a large number of wind 
farms or PV systems, if the demand response is applied using 
an appropriate price structure. 

I. Introduction 

To mitigate the effect of climate change, it is essential 
to reduce the CO2 emissions in various industry sectors 
by using renewable energies. The deployment of renewable 
energies such as wind and photovoltaic (PV) systems has 
begun in the power industry sector. The share of renewable 
energies in the total electricity generation will reach 15% 
in the near future. For instance, the climate change package 
"triple 20" is scheduled to take effect in the European Union; 
the implementation of this package is expected to reduce 
greenhouse gas emissions by 20% by the year 2020 compared 
to 1990 levels, cutting consumer energy consumption by 20% 
by improving energy efficiency, and ensuring that 20% of the 
EU's energy mix comes from renewables. 

It is however noted that the output of renewable energies 
fluctuates significantly depending on weather conditions. With 
an increasing share of energy coming from the renewable 
energy sources the share of energy contributed by various 
types of thermal power plants is expected to decrease. At the 



same time, there is an increasing fluctuation in the supply of 
electric power. This implies that the supply-demand balance 
of electric power cannot be attained using the currently used 
conventional power systems. 

Smart grids have the potential to increase the capability 
to attain the supply-demand balance of electric power by 
using new technologies, such as smart meters, home energy 
management systems, building and energy management sys- 
tems, control of charge and discharge of electric vehicles, 
and various types of electric storages |T|, f2l, lO, A 
power utility develops a plan for the operation of a power 
system based on the forecast output of renewable energies 
for 24 h. The forecast output includes time series of output 
scenarios and its error for the aggregate outputs of wind power 
and PV system in the power system under consideration. 
An accurate forecast makes it possible to operate the power 
system economically without making an overestimation of the 
required adjustments to meet the demand and supply. 

From our studies on the operation cost, the CO2 emissions 
and the peak load of the smart grid fSl, ||6|, we identified 
the need to have a detailed understanding of the time evo- 
lution of the power system. In this paper, we describe a 
unit commitment model [7|,[8|,|9|,|10|,[11|,|12| with demand 
response for the integration of renewable energies to analyze 
requirements of the forecast output and its error. The paper 
is organized as follows. In Section II, the formulation of 
the model is explained. In Section III, the analysis of a 
power system is described. Finally, Section IV presents our 
conclusions. 

II. The Model 

The purpose of our unit commitment model is to obtain the 
time series for the operational state of thermal power plants 
that would maximize the profit of an electric power utility by 
taking into account both the forecast of output and its error for 
renewable energies |lT3l, Oil, |[l3, |[l6l, |[l7l and the demand 
response of consumers on the change of electricity prices 1 18|, 
OSl, El, im, lEl, El. The model is formulated as a mixed 
integer linear programming problem. 

A. Objective Function 

The time series of the operational state of thermal power 
plant i{i = 1, • • • , N) is obtained by maximizing the objective 
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Fig. 1. Demand and supply 



function: 
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This objective function represents the profit of an electric 
power utility. The first term of the r.h.s. in Eq. ([T]) is the 
sales revenue and the second term is the operation cost. 
Here A'', T, L are the number of thermal power plants, 
time horizon, and the number of price levels, respectively. 
Continuous variables p\ is the output power variable of thermal 
power plant i, and integer variables u\, zj, w\ are the status 
production variable of thermal power plant i (l=committed, 
O=decommitted), the start-up variable of thermal power plant 
i (l=start up, O=others), and the demand response variable 
(l=selected, O=not selected), respectively. Parameters Si and 
fo^are the start-up cost of thermal power plant i and the fuel 
cost of the thermal power plant i, respectively. The forecasted 
demand and its error are indicated by dp' and a^, respectively. 
Here (/) stands forforecasting. 

Other parameters related to the demand response f, r\ 
are the average electricity price, the price level, and the price 
elasticity of demand, respectively. If the electricity price r 
deviates from the average price r, the demand d is changed 
from the average demand d as follows: 
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The dependence of the demand d on the price r is depicted 

in Fig. m i24i, ca, Ea, eu, ed. 

B. Global Constraints 

The sum of the demand response variable wl has to satisfy 
the constraint 
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to ensure that only a single price level r' is selected. In addition 
to this constraint, the average of the selected price r' has to 



be equal to the average price r 

T L 



(4) 
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Moreover the total demand has to be unchanged by the demand 
response: 

yd[fUTd\f^yyJir\ (5) 
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Importantly, the sum of supply has to be greater than the 
demand: 



N 
i=l 



wd, 



if) 



-pv[f~>>dY'Y.^\ 

1=1 



(6) 



where wdY ' and pv'"/^ are the forecasted wind power gener- 
ation and the forecasted PV generation, respectively. 

If we consider the forecast eiTor of demand ad, the forecast 
eiTor of wind power cr^, and the forecast error of PV cTp, the 
constraint in Eq.® can be rewritten as 

T.tiP\ + ^d[f^ +P^t'' 



(7) 

where a and are the probability to ensure the supply- 
demand balance and the cumulative distribution function, 
respectively. If the forecast error is distributed according to 
the normal distribution, the probability density function is 
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and the cumulative distribution function is written using the 
error function erf[-] as 



x) = - [1 + erf 



(9) 



where ji and cr are a mean and a standard deviation, respec- 
tively. However, if the probability density function p{x) is a 
Laplace distribution 
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then the cumulative distribution function 

2' 

Here a standard deviation is given by cr = \/2b and sgn{x — 
n) = -\-{x > fi),—{x < /x). The functional forms for these 
distributions are depicted for p, — and cr = 1 in Fig. |2] It is 
noted that here the Laplace distribution shows a distribution 
tail longer than the normal distribution. This implies that the 
probability to ensure supply-demand balance differs in these 
two distributions. 

C. Local Constraints for Thermal Power Plants 

The following constraints are used for each thermal power 
plant as typical constraints in a unit commitment model. 
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Fig. 2. Distribution of forecast error 



1) Generation Capacity: The output power has to be 
between the maximum output power p^-^g^^ and the minimum 
output power p^j„ when the operation is in steady state: 



(12) 



2) Ramp-up Limit: The increase in the output of thermal 
power plant i has to be smaller than the maximum ramp-up 
speed A_|_ when the unit is up at the previous time step and is 
smaller than the minimum output power p^j„ when the unit 
is down at the previous time step: 



(13) 



3) Ramp-down Limit: The decrease in the output of thermal 
power plant i has to be smaller than the maximum ramp-down 
speed A_ when the unit is up at time step t and is smaller 



than the maximum output power 
at time step t: 



when the unit is down 



(14) 



4) Minimum Up-time Constraint: Thermal power plant i 
has to be operated longer than the minimum up-time require- 
ment T^, once the unit is up: 
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5) Minimum Down-time Constraint: Thermal power plant 
i has to be stopped longer than the minimum down-time 
requirement ri, once the unit is down: 
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TABLE I 

The parameters of thermal power plants 



i 


bi 


Pmax 


p' - 


A+ 


A_ 




ri 


St 


1 


3.0 


50.0 


25.0 


0.5 


0.5 


3.0 


3.0 


1000.0 


2 


3.0 


50.0 


25.0 


0.5 


0.5 


3.0 


3.0 


1000.0 


3 


3.3 


15.0 


7.5 


0.5 


0.5 


3.0 


3.0 


200.0 


4 


3.3 


15.0 


7.5 


0.5 


0.5 


3.0 


3.0 


200.0 


5 


4.3 


6.0 


2.0 


5.0 


5.0 


3.0 


3.0 


100.0 


6 


4.3 


6.0 


2.0 


5.0 


5.0 


3.0 


3.0 


100.0 


7 


4.3 


6.0 


2.0 


5.0 


5.0 


3.0 


3.0 


100.0 


8 


7.1 


10.0 


4.0 


5.0 


5.0 


3.0 


3.0 


200.0 


9 


7.1 


10.0 


4.0 


5.0 


5.0 


3.0 


3.0 


200.0 


10 


7.1 


10.0 


4.0 


5.0 


5.0 


3.0 


3.0 


200.0 


11 


9.0 


5.0 


2.5 


0.5 


0.5 


3.0 


3.0 


100.0 


12 


9.0 


5.0 


2.5 


0.5 


0.5 


3.0 


3.0 


100.0 



6) Constraint on the Start-up Variable: The start-up vari- 
able zl has to satisfy the following constraints by definition: 



uU{t > 2). 
III. Analyzed Power System 
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We analyze the time series for the operational state of 
thermal power plants by taking into account both forecast 
output and its error for renewable energies and the demand 
response of consumers. We consider a small power system 
consisting of thermal power plants, PV systems, and wind 
farms in this study, although the model is easily extendible to 
a larger power system. The number of thermal power plants is 
12, and parameters are given in Table U The installed capacity 
for each PV system and wind farm is 30MW. The intra-day 
peak demand is about 170MW. We assume the forecast errors 
cr^ and (Tp are both 10%(3MM^) of the installed capacity 
in the reference case. In addition to the reference case, we 
analyze two more cases a.w = 6MW and cr^ = 9MW, while 
(Tp remains unchanged. The forecast error for the demand is 
ad — for all cases. Scenarios for the demand, the PV output 
power, and the wind output power in the reference case are 
shown in Fig. [3] While wind scenario 1 is used in the reference 
case, wind scenario 2, where wind power suddenly decreases 
in the evening, is used to discuss the fluctuation of wind output 
and the effect on the demand response. 

The marginal cost C^^ , the spinning reserve i?f , and the 
load frequency control (LFC) margin i?f are calculated using 

, N). The marginal cost 



the obtained time series pl{i = 1. 
Ct" is defined by 

= max{5i|uj = 1} 



(18) 



The spinning reserve i?f and the LFC margin are defined 
as proxy quantities in this study as follows: 
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Fig. 3. Given load, PV, wind scenarios 
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Here, p). is the output power of thermal power plant i obtained 
in the optimization with the constraint of Eq.©. The spinning 
reserve i?f is ready to generate power to absorb the fluctuation 
instantaneously. On the other hand, the LFC margin is the 
remaining capacity that is able to increase output with the 
constraints of Eq. ([13]) and Eq. (fT4b . 

IV. Results 

We discuss the requirements on the forecast output and 
its error for the integration of the renewable energies in this 
section. The unit commitment model formulated as a mixed 
integer linear programming problem in Eqs. ([T]i to ( fTTI ) was 
solved to analyze the power system described in Section [III] 
using a commercial solver ll29l . If30l . 

A. Reference Case 

The results of the reference case are shown in Fig.lD Figure 
m (a) depicts the share of the thermal power plants, the wind 
power, and the PV systems to satisfy the given demand. The 
thermal power plants are numbered in the order of increasing 
operation cost, i.e., the merit order. The fluctuation of the 
wind and the PV outputs was absorbed by starting the thermal 
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Fig. 4. Results of the reference case 



power plants serially in the merit order. Figure |4] (b) shows 
that the mai-ginal cost C^' is high from 9 : 00 to 20 : 00. 
This is because some additional thermal power plants are 
operating in this period. The price structure where the daytime 
rate is higher than the nighttime rate is obtained as a result 
of maximizing the profits of an electric power utility. In a 
reflection of the price structure, the sales revenue is high 
during daytime as shown in Fig.|4](c). Figure |4](d) depicts that 
the spinning reserve i?f has a flat structure during daytime, 
while the LFC margin R^ is high during nighttime. 

B. Demand Response 

It is well known that the price elasticity of demand 
is small because the share of expenditure for electricity in 
a household budget is low ISTl . The effects of the demand 
response with = —0.30 are shown in Fig. |5] Figure |5] (a) 
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Fig. 5. Demand response 



depicts that the load profile becomes flat and the peak load 
is reduced by about 2QMW. The profile of the marginal cost 
becomes flat and the price structure is significantly changed as 
shown in Fig. |5] (b). In general, a larger effect to the demand 
response is expected for a larger e^. 

C. Effects of the Forecast Error 

The effects of the forecast error on the operation cost and 
reserves are analyzed here. The forecast errors and Cp are 
\Q%{iMW) of the installed capacity in the reference case. In 
addition to the reference case we analyzed two more cases; 
f u) = CTp = QMW and = cFp — 9MW. Figure |6] (a) 
depicts that the operation cost increases with the forecast error 
(Ttu. This tendency is also true for the spinning reserve i?f 
as shown in Fig. |6] (b). The intra-day structure is the same 
for the three cases, but the level increases by 5AMW as 
(Ttu increases by 3MW. This is reasonable because a larger 
reserve is required to absorb the larger fluctuation. Figure |6] 
(c) shows that the LFC margin i?f deer eases during daytime, 
while the margin for all three were relatively high during 
nighttime, independently of the forecast error This implies that 
the demand during nightime is relatively low, and therefore, 
the remaining capacity is large in the operating thermal power 
plants. 

D. Effects of Sudden Decrease in Wind Power 

In wind scenario 2, the wind power suddenly decreased in 
the evening compared to wind scenario 1 . We expect that more 
thermal power plants will be operated to generate power to 
absorb this sudden decrease in wind power The effects of the 
sudden decrease in wind power were analyzed and the results 
are shown in Fig. [7] Figure [T] (a) depicts an increase in the 
number of operating thermal power plants between 17 : 00 
and 20 : 00. This increase is more clearly shown in Fig. |7] 
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Fig. 6. Effect of forecast en'or on cost and reserves 



(b). Nine thermal power plants are in operation during 17 : 00 
and 20 : 00 for scenario 2, whereas just seven plants operate 
in the same period for the scenario 1. The increase in the 
number of operating thermal power plants in such a short 
period does not affect the total operation cost significantly. 
It is however noted that the power utility has to continue to 
use thermal power plants for ensuring supply-demand balance 
even after installing a large number of wind farms or PV 
systems. Therefore, the substitution of thermal power plants 
by wind farms or PV systems is not expected to be very high, 
although this issue has to be studied quantitatively using actual 
data. 

E. Demand Response and a Sudden Decrease in Wind Power 

Finally, the effects of the demand response with = —0.30 
in wind scenario 2 were analyzed. If the number of operating 
thermal power plants is reduced by the demand response, the 
economic value of the wind farms will increase due to the 
high number of thermal power plants substituted by wind 
farms. Figure [8] (a) shows that the number of operating thermal 
power plants is reduced by the demand response during 9 : 00 
and 20 : 00. This implies that the effects of the demand 
response can be clearly observed in the number of operating 
thermal plants. Figure [8] (b) depicts the price structure with 
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Fig. 7. Wind scenario 2 



and without the demand response. The rectangular structure 
for wind scenario 2 is smoothened by the demand response. 
The price during 17 : 00 and 20 : 00 is kept at a high level. 
It is recognized that the marginal cost is high during 17 : 00 
and 20 : 00 in Fig. [8] (c). Therefore, the high price during 
17 : 00 and 20 : 00 is a consequence of the high marginal 
cost during the same period. Because the operation cost was 
not increased in scenario 2, we cannot expect the demand 
response to cause a reduction in the operating cost in this 
case. However, a power utility has to continue to use thermal 
power plants for ensuring supply-demand balance; some of 
these plants can be decommissioned after installing a large 
number of wind farms or PV system, if the demand response 
is applied using an appropriate price structure. 

V. Conclusion 

The output of renewable energy fluctuates significantly 
depending on weather conditions. Therefore, it will be difficult 
to ensure supply-demand balance of electric power using cur- 
rently used conventional power systems. We have developed 
a unit commitment model to analyze requirements of the 
forecast output and its error for renewable energies. Our model 
obtains the time series for the operational state of thermal 
power plants that would maximize the profits of an electric 
power utility by taking into account both the forecast of output 
and its error for renewable energies and the demand response 
of consumers. The model is formulated as a mixed integer 
linear programming problem. 

We considered a power system consisting of thermal power 
plants, PV systems, and wind farms. The unit commitment 
model was solved using a commercial solver to analyze the 
power system. 

First, the basic property of the model was discussed using 
the results of the reference case. The fluctuation of the wind 
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Fig. 8. Wind scenario 2 with demand response 



and the PV outputs was absorbed by starting the thermal power 
plants serially in the merit order. The price structure where the 
daytime rate was higher than the nighttime rate was obtained 
as a result of maximizing the profits of an electric power utility. 

Next, the effects of the forecast error on the operation cost 
and reserves were analyzed. In addition to the reference case 
((Tto = 3MW), we analyzed two more cases, namely, a^j ~ 
(jp = 6MW and aw — (Jp = 9MW. We confirmed that the 
operation cost increased with the forecast error The intra-day 
structure of the spinning reserve was the same for the three 
cases, but the level increased with the forecast error. The LFC 
margin decreased during daytime, while it remained relatively 
high during nighttime, independently of the forecast error 

Then, the effects of a sudden decrease in wind power were 
analyzed. More thermal power plants will have to be operated 
to generate power to absorb this sudden decrease in wind 
power The increase in the number of operating thermal power 
plants within a short period did not affect the total operation 
cost significantly; however the substitution of thermal power 
plants by wind farms or PV systems is not expected to be very 
high. 

Finally, the effects of the demand response in the case of a 
sudden decrease in wind power were analyzed. We confirmed 
that the number of operating thermal power plants reduced 
by the demand response. This implies that the number of 



operating thermal power plants is controlled efficiently by the 
demand response if an appropriate price structure is used. A 
power utility has to continue to use thermal power plants for 
ensuring supply-demand balance; some of these plants can be 
decommissioned after installing a large number of wind farms 
or PV systems, if the demand response is applied using an 
appropriate price structure. 

In future work, we intend to study the demand response 
technology and the forecast output of wind farms and PV 
systems using the developed unit commitment model for the 
maximizing the renewable energy integration in an actual 
power system. We also plan to quantitatively study the issue 
of substitution of existing thermal power plants by renewable 
energy using a large set of relevant data ll32ll . 
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